ID-R148  396 SOLI TONS  AND  SEASAKU)  NITRE  CORP  HCLEflN  VA  K  I  CASE  1/1 
AUG  84  JSR-84-202D  F19628-84-C-0001 


UNCLASSIFIED 


F/G  20/4 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS  ~  1965  “  A 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Oau  Entaradl 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

JSR-84-203D 


4.  TITLE  (and  Subtitla) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


SoliCons  and  SeaSat 


6.  PERFORMING  ORG.  REPORT  NUMBER 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

The  MITRE  Corporation 
1820  Dolley  Madison  Blvd. 

McLean,  VA  22102 


11.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 


8.  CONTRACT  OR  GRANT  NUM8ERU) 


F19628-84-C-0001 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE  13.  NO.  OF  PAGES 

August  1984  18 


15.  SECURITY  CLASS.  <o«  thi«  report) 


14.  MONITORING  AGENCY  NAME  &  ADDRESS  (if  diff.  from  Controlling  Office) 


Unclassified 


15e.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thi*  report) 


distribution  statement  a  j 

Approved  lot  pubHc  nImmi  I 
DMiMiot  Unlimited  I 


17.  DISTRIBUTION  STATEMENT  (of  th«  abstract  antaracf  In  Block  20.  If  difforont  from  raport) 


19.  KEY  WORDS  (Continue  on  ravarta  sida  if  nacatsary  and  identify  by  block  number) 


20^  ABSTRACT  (Continue  on  reverse  sida  if  necessary  and  identify  by  block  number) 

It  has  been  suggested  that  Soliton  formation  might  be  relevant  to  SeaSat  obser¬ 
vations.  It  has  also  been  said  that  there  are  no  Solitons  in  more  than  one 
space  dimension.  The  second  statement  is  demonstrated  to  be  false.  The 
Kadomtsev- ’etviashvile  equation  relevant  to  Internal  Waves  is  shown  not  to 
have  SoliL  .  solutions.  This  lends  support  to  the  view  that  Solitons  and 
SeaSat  have  little  in  common. 


DD°ri473 

EDITION  OF  1  NOV  65  IS  OBSOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  this  PAGE  ivNPe"  Data  Emeradi 


Solitons  and  SEASAT 


K.  M.  Case 


August  1984 
JSR-84-203D 


Approved  for  public  release,  distribution  unlimited 


JASON 

The  MITRE  Corporation 
1820  Dolley  Madison  Boulevard 
McLean.  Virginia  22102 


DTIC 

SFLECTE 
DEC  1  0  D84 

B 


ABSTRACT 


I 

It  has  been  suggested  that  Soliton  formation  might  be 
relevant  to  SeaSat  observations.  It  has  also  been  said  that  there 
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I  are  no  Solitons  in  more  than  one  space  dimension.  The  second 

statement  is  demonstrated  to  be  false.  The  Kadomtsev-Petviashvile 
equation  relevant  to  Internal  Waves  is  shown  not  to  have  Soliton 
solutions.  This  lends  support  to  the  view  that  Solitons  and  SeaSat 
have  little  in  common. 
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1.0  INTRODUCTION 


The  persistence  of  long  V-shaped  wakes  observed  by  SeaSat 
has  led  to  suggestions  that  the  phenomena  may  be  related  to  Internal 
Wave  Solitons.  Most  observations  were  made  under  conditions  for 
which  one  would  have  little  or  no  reason  to  expect  Solitons  to  be 
relevant.  However,  there  is  one  case  (Rev.  407)  for  which  one  might 
think  otherwise.  This  was  a  big  ship  in  shallow  water  with  a  strong 
thermodine. 

A  priori  there  is  little  likelihood  that  the  equations 
describing  internal  waves  generated  by  ships  will  admit  Soliton 
solutions.  (By  a  Soliton  we  mean  a  non-singular  disturbance 
localized  in  space  at  any  time  which  retains  its  integrity  on 
interactions  with  similar  disturbances.)  While  there  are  many 
equations  in  one  space  dimension  which  have  Soliton  solutions,  little 
is  known  about  spaces  of  higher  dimensions.  (Indeed  it  is  sometimes 
said  that  there  are  no  true  solitons  in  case  the  number  of  space 
dimensions  is  greater  than  one.  Below  we  will  see  this  is  not  true.) 

Since  the  paradigm  of  an  equation  describing  Solitons  is  the 
Korteineg-deVr ies  (KdV)  equation  it  is  natural  to  take  as  our 
starting  point  an  equation  closely  related  to  this  which  does 


describe  fully  three  dimensional  internal  waves.  (The  class  of 


solutions  described  is  not  exactly  those  we  would  expect  to  be  ship 
produced.  However,  it  is  hoped  that  the  results  obtained  will  give 
support  to  our  conviction  that  Internal  Wave  Solitons  are  not 
relevant  for  the  SeaSat  photographs.) 


where 


The  equation  we  have  in  mind  is 


u  +  3  *  u  +3  (3u^  +  u  ) 

t  x  yy  x  xx 


X  00 

3’S  -  V2  {  /  *  /  }  ♦  <*'>  dx' 

x 


Closely  related  to  Eq.  (1)  is  the  equation 


u  -  3  *  u  +3  (3u^  +  u  )  ■  0 

t  x  yy  x  xx 


The  Eqs.  (1)  and  (2)  are  known  as  the  Kadomtsev  - 
Petviashvili  (l)  equations.  A  derivation  of  Eq.  (1)  in  the  case  of 
internal  waves  is  given  in  reference  (2).  An  important  point  is  that 
for  Internal  waves  the  sign  of  the  term  in  Eq.  (1)  involving  Uyy  is 
unambiguously  required  to  be  plus.  (There  are  physical  situations 


2 


where  Eq.  (2)  can  hold.  An  example  ia  when  capillarity  is 
important.) 

Our  main  result  is  the  following:  The  extension  of  the 
N-Soliton  solutions  of  K-deV  to  solutions  of  Eqs.  (1)  and  (2) 
generally  break  up  asymptotically  into  K-deV  Solitons  moving  in 
arbitrary  directions  in  the  x-y  plane.  (These  are  then  plane 
waves.)  They  are  physically  unacceptable  as  being  non-localized. 
However,  for  special  values  of  the  parameters  these  solutions  are 
localized  in  x,y.  However,  for  Eq.  (1)  these  localized  solutions  are 
singular.  Hence  these  solutions  again  are  not  physical.  (This  is 
not  so  for  Eq.  (2).  True  localized,  non-singular,  non-interacting 
"lumps"  result.  We  include  these  primarily  to  show  that  multi¬ 
dimensional  Solitons  do  indeed  exist.) 

In  Section  2  the  formalism  used  by  Zakharov  and  Shabat  (3) 
to  integrate  the  K-P  equations  is  summarized.  The  generalization  of 
the  N-Soliton  KdeV  solution  Is  given  in  Section  3.  In  Section  4  the 
"lumps"  which  result  when  special  relations  exist  between  the  N- 
Soliton  parameters  is  presented. 


”  • 

2.0  FORMALISM 

”  • 

Zakharov  and  Shabat  (3)  have  introduced  a  method 

to 

'  .'  *  __  " 

integrate  Eqs.  (1)  and  (2).  We  summarize  this  here. 

• 

Let  F(x,z;y,t)  satisfy  the  two  equations 

-  »■ 

4  3t  3xJ  3zJ 

(3) 

• 

and 

~  • 

2  2 

I  3F  .  3  F  3  F  . 

(4) 

• 

Determine  K(x,z;y,t)  from  the  Volterra  equation 

00 

F(x,z)  +  K(xfz)  +  /  K(x,s)  F(a,z)ds  *  0 

X 

(5) 

(Here  we  have  suppressed  the  parametric  arguments  y  and  t.) 


Then  they  show  that 


u  -  2  —  K(x,x) 

dx 

satisfies  Gq.  (1). 


(6) 


To  find  solutions  of  Eq.  (2)  we  merely  note  that  from  any 
solution  of  Eq.  (1)  we  can  obtain  a  solution  of  Eq.  (2)  by  the 
replacement  y  ♦  iy  . 
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3.0  GENERALIZATION  OF  THE  N-SOLITON  K-deV  SOLUTIONS 


We  follow  reference  (3). 


A.  Suppose 


F 


-KTX 

e 


nz 


M(t,y) 


From  Eq. 


M 


Then  Eq. 


C 


(3)  we  find 

N  4(k3  +  n3)  t 
*  C(y)  e  ' 


(4)  gives 

-  M  e_/J  (<2  "  n2)  y 
o 


* 


* 


i.e 

[/3  (r,2  -  <2)y  +  4(k3  +  n3)t] 

M  -  M  e 
o 


Our  Gelfand-Levitan  Equation  (5)  becomes 


Clearly  K(x,z)  -  K(x)  e  1 


Inserting  in  Eq.  (11)  yields 


i//  \  u  -<x  ,  H  e 
K(x)  +  Me  +  - 


■(<  +  n)x 


K  =  0 


(12) 


<  +  n 


Solving  gives 


K(x,z) 


-M  e 


-(<x  +  nz) 


(13 


,  ,  M  -(<  +  n)x 
1  +  e 


Then 


K(x,x) 


-M  r(<  *  n)X 

1  +  -5-  e’(‘  +  ")x 
<  +  n 


^  u  (i  + 


-(<  +n)x 


«  +  n 


and  so 


£ 

3x* 


2  - —  In  [l  +  -  ■ 

2  <  +  n 


-(<  +  n)x 


Defining  xQ  by 


(14 


gives 


M 


In 


+  /3  (n  -  <)y  +  4(<  -  <n  +  n  )t  (15) 


k  +  n  k  +  n 


Then  Eq.  (14)  becomes 


V?  (<  +  njL 


cosh2  (x  “  V 


(16) 


To  interprete  this  consider  k  =  n  then 


2< 


cosh  ic(x  -  x  ) 
o 


(17) 


with 


x  =  constant  +  4k  t 
o 

This  is  just  a  K-deV  Soliton.  For  k  ^  n  the  solution  (16) 
is  then  a  K-deV  type  soliton  propagating  at  an  arbitrary 
angle  with  respect  to  the  x-axis.  This  is  a  plane  wave  -  it 
is  constant  on  the  line 

2  2 

x  -  /3  (n  -  k  )y  *  constant  , 
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and  hence  physically  really  not  acceptable. 

The  above  is  readily  extended  to  get  the  2-dimensional 
extension  of  the  N-Soliton  K-deV  solution.  Thus  we 
generalize  Eq.  (10)  by  choosing 


-<  x  -  n  z 

F  =  Z  Mn  (t ,y)e  n  n 
n 


(ii 


Since  Eq.  (3)  and  (4)  are  linear  we  have  in  analogy  to  Eq 

(10) 


M  «*  M  (o)  e 
n  n 


[/3  (tt  -  tc^)y  +  4(x^  +  )t] 

n  n  n  n 


(1' 


The  degenerate  integral  equation  for  K(x,z)  is  then 
satisfied  by 

-n  z 

K(x,z)  *  Z  Kn(x)  e  "  (2( 


where 


-<  x  -( <  +  n  )x 

K  (x)  +  M  e  n  +  M  Z  - - - - —  K  (x)  »  0  ( 2 

n  n  n  ,  mv 

m  k  +  n 


n 


m 


This  is  readily  solved  using  Cramer's  rule.  Thus 


Kn(x) 


Vx) 

A(x) 


(22) 


where 


A  =  det  {6  +  — 2- 

1  nm 


m  +  n  )x 

Men  m 


(23) 


<  +  n 
n  m 


and  Ar  is  obtained  from  A  by  replacing  the  n'th  column  of 
the  matrix  by  the  vector 


-K  X  -K?X 

(-M^e  ,  -M2e  ,  ...) 


Suggested  by  the  form  of  Eq.  (14)  we  look  at  E  A  (x)  e 

n 

and  verify  that 


L  An(x)  e  nX  =  A(x) 
n 


-nx 


therefore  K(x,x)  *  -  In  A  and  so 

9x 


u(x)  *  2  -  In  A 

3x2 


(24) 


To  interpret  this  result  we  consider  the  limit  as  y,  t 
go  to  infinity.  Assuming  no  special  relations  between  the 
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various  pairs  (k  ,n  )  .  Then  we  note  chat  for  large  |y|, 
n  n 

|t|  there  will  be  regions  where  one  of  the  Mn  (say  Mj)  is 
much  larger  than  all  others.  Hence  then 


H,  .-(‘j  +  V* 
A  •  1  +  — J - 

V*J 


(25] 


which  is  just  the  single  soliton  result  of  Eq.  (14).  Thui 
asymptotically  the  solution  breaks  up  into  a  sum  of  the 


4.0  DEGENERATE  CASES 


It  is  well  known  that  the  K-deV  equation  has  in  addition  to 
solutions  like  those  of  Eq.  (24)  solutions  which  are  rational 
functions  of  the  coordinates.  These  can  be  obtained  by  making  a 
suitable  ansatz  for  the  form  of  solution.  This  is 


U  *  £  - - -  (26) 

11  [x  -  an(t)f 

Inserting  in  the  K-deV  shows  that  this  will  be  a  solution  if  a  ■  2  , 
the  cR  are  constants  and  the  aR  satisfy  simple  coupled  ordinary 
differential  equations.  An  alternative  approach  is  to  take  the 
general  solution  of  Eq.  (24),  specify  relations  between  the 
parameters,  *n,  nn»  and  pass  to  limits.  This  is  the  procedure  given 
in  reference  (4)  -  and  the  one  we  will  follow. 


We  have  seen  that  a  solution  is  obtained  from  Eq.  (24) 
with  A  given  by  Eq.  (23).  Introduce  X^,  Yn  by 


X  +  Y 
n  n 


X  -  Y 
n  n 


then  choosing  M  (0)  ■  -a  X  we  have  A  ■  det  T  where 
n  n  n 
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-fx  +  y  +  X  -  y  lx 
1  n  n  m  mJ 


r  -  5  + 

nn  nn 


Y  -  Y  +  X  +  X 
n  m  n  m 


X  [  ] 
nl  Jn 


M  »  -  a  X  e 
n  n  n 


M„  -  h*  v  +  (A*  ♦  3^,t] 


Now  look  at  the  limit  as  all  X  +0 

n 


Assuming  a  1  +  5  X  we  have  on  expanding  In  X  and 
n  n  n  n 

keeping  only  matrix  elements  of  first  order  In  the  X 


r  -  x  d 

nn  n  n 


n  *  m  T  ■ 
nm 


Y  ~  Y 
n  m 


where 


d  -  x  -  C  +  /3  y  y  '  3y  t 
n  n  n  n 


.N.  , 


det  r  -  (nx  )  det  r 
n 


with  T'  -  d 
nn  n 


n  *  m  V 


We  can  then  write 


y  -  y 

n  m 


u  -  2  — —  in  det  T' 

.  2 


The  essential  points  are  seen  by  considering  the  case  of 


N  »  2.  Then 


A'  -  det  T'  ■  d^  dj  + 


<T1  -  V 


*  [x  -  -  3Yjt][x  -  «2  +  /lY2y  -  3Y2t] 


(Y,  -  Y2) 

Clearly  u  is  then  a  rational  function  of  x,  y,  t.  Consider  various 

possibilities  for  C  ,  y  . 

n  n 


(i)  C  ,  y  real.  The  solution  is  clearly  singular, 
n  n 


1A 


(There  are  points  where  A'  •  0  .) 


(ii)  c  ,  y  complex.  For  the  solution  to  be  real  we  must 
n _ n _ 

have 


* 


Then 


A' 


lx 


c  +  /3  Yy  -  3y2t 1 2 


(38) 


The  second  term  can  vary  from  0  to  "  and  therefore  A' 
vanishes  at  least  once. 


These  rational  solutions  of  Eq.  (1)  are  physically 
unacceptable. 

On  the  other  hand  we  noted  that  a  solution  of  Eq.  (2)  is 
obtained  from  one  of  Eq.  (1)  by  the  replacement  y  +  iy  .  In  such  a 
case 


A'  ♦  [x  -  c1  +  i/3yxy  -  3y*t][x  -  C2  +  i^2y  "  3y2c^ 
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5.0  CONCLUSIONS 


Two  dimensional  Solltons  exist.  However,  some  equations 
describing  Internal  waves  which  might  be  strongly  suspected  of  having 
Soliton  solutions  do  not. 

While  not  definitive  this  lends  significant  support  to  the 
view  that  Solltons  have  no  connection  with  the  phenomena  observed  by 


SeaSat. 
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